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A method is proposed to test for the nature of the pseudogap phase in cuprates using the recently 
developed technique of Fourier transform scanning tunneling spectroscopy. We show that the ob- 
served quasiparticle interference patterns depend critically on the quasiparticle coherence factors, 
making it possible to distinguish between the pseudogap dominated by superconducting fluctuations 
and by various particle-hole condensates. 



Existence of the pseudogap phase, a non- 
superconducting phase with suppressed single particle 
density of states (DOS) represents perhaps the most 
dramatic departure of cuprates from the Landau Fermi 
liquid BCS-Eliashberg paradigm believed to describe 
conventional low-T c superconductors. Theoretical 
approaches to the pseudogap phenomenon in cuprates 
fall broadly into two classes. One school of though 
attributes it to the incipient superconducting order 
whose amplitude forms at T* > Tn but remains phase 
incoherent down to T c @, S, 0, H II 0, S 0- The 
other school ascribes the pseudogap to the formation 
of some other (static or fluctuating) order, usually in 
the particle- hole (p-h) channel EH El El El • Progress 
in understanding the physics of cuprates, and perhaps 
other strongly correlated superconductors, depends 
on the successful determination of the origin of the 
pseudogap phenomenon. 

In this Communication we propose a test for the na- 
ture of the pseudogap phase based on the study of the 
quasiparticle interference patterns that can be observed 
by means of Fourier transform scanning tunneling spec- 
troscopy (FT-STS) El El [l6|. The existing data, taken 
deep in the superconducting phase of E^S^CaC^Os+tS 
(Bi-2212), exhibit distinctive patterns of peaks in the re- 
ciprocal space. These peaks disperse through the Bril- 
louin zone in a manner consistent with the band struc- 
ture deduced from the angle resolved photoemission spec- 
troscopy (ARPES). In what follows we demonstrate that, 
contrary to the existing consensus, the patterns observed 
in FT-STS depend in a crucial way on the quasiparticle 
coherence factors. We show that the appearance of peaks 
at the particular positions is unique to the superconduct- 
ing order, and reflects the pairing (p-p) character of the 
coherence factors. By contrast the interference patterns 
in a state with p-h ordering will be qualitatively different 
in that peaks appear at different k-vectors or not at all. 
Our proposal consists of extending the FT-STS measure- 
ments into the pseudogap phase above T c . Observation 
of patterns similar to those seen well below T c would then 
imply that pseudogap is predominantly of superconduct- 
ing origin. Observation of qualitatively different patterns 
discussed below would imply order of another type. Pre- 
liminary data E3 appear to support the former scenario. 

The existing FT-STS results El E3 have been inter- 
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FIG. 1: Right: The banana shape contours of constant low 
energy around the nodes of the gap and seven vectors connect- 
ing their ends [141] . Left: Numerical evaluation of the Fourier 
transformed local density of states in Eq. (J with t' — 0.3t, 
Ao = 0.2t, fi = — t and u — 0.45Ao. In order to simulate 
the experimental resolution we filtered out the low-intensity 
background signal. The peaks disperse in a manner consistent 
with the octet model and the experimental data [lH IT3|. 



preted via the 'octet model' illustrated in Fig. [I] 
This model asserts that the interference patterns arise 
due to the elastic quasiparticle scattering from random 
disorder between the regions in the Brillouin zone with 
high DOS. In a d-wave superconductor (dSC) these are 
situated at the ends of the banana-shaped contours of 
constant energy and lead to seven characteristic vectors 
q^. This octet model works remarkably well in describing 
the data and furthermore agrees with detailed numerical 
studies of the interference patterns 0, 0, [2(J . Our first 
step will therefore be to understand how the interference 
pattern is formed in the superconducting phase. 

The local density of states (LDOS) in a superconductor 
n(r,w) is given by the particle-hole part of the Nambu 
Green's function. Within the usual T— matrix formula- 
tion |l8|, LUj the Fourier transformed LDOS modulation 
5n(q_,u) can be expressed as 



Sn 



lm j2 ^[Go(k,o;)f(k,k - q,uOG (k - q,o;)]ii, 



{1) 

where L is the linear size of the system and Go(k, iuo) = 
[iuo — 6kT3 — AkTi] -1 is the superconducting Green's func- 
tion in the Nambu space. In what follows we consider 
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a model with dispersion ek = —2t(cosk x + cosk y ) — 
4t f cos k x cos k y — \i and a d-wave gap Aw = ^ (cos k x — 
cos k y ). As shown by Capriotti et al [2(^ the interfer- 
ence pattern measured in FT-STS for weak non-magnetic 
impurities is encoded in the quantity 



A(q,cj) = ^[G (k,cj)r 3 Go(k-q,cj)]ii. 



(2) 



For simplicity we start by analyzing this case and we 
return to Eq. Q with the full T-matrix shortly. 

It is instructive to evaluate A(q, uj) in the low energy 
limit a; <C A , where analytic results can be obtained. 
Inserting Go into Eq. (0) we find 

1 (iw + + e_) — A+A_ . N 



with e± = e k±q/2 , A± = A k±q/2 and E± = ^ e 2 ± + A 2 ± . 

We now focus on the situation when both k± |q are close 
to one nodal point; this should contribute to the peak in 
A(q, uj) at the vector labeled as qj in the Fig. [l] Near 
this node we define a local coordinate system (fci,fc2) 
and linearize the dispersion in the usual manner; ek — > 
VFki and Ak — > ^a&2, where and ^a are quasiparticle 
velocities perpendicular and parallel to the Fermi surface, 
respectively. We find 
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(k-q) 2 ][^ 2 + (k + q) 2 ]' 
(4) 

where we have scaled the integration variables k\ — > 
ki/vp, k2 — > fe/^A and defined a new vector q = 
^(v F q 1 ,v A q2)- 

The above rescaling leads to our first useful insight: in 
the scaled frame of reference the contours of constant en- 
ergy are concentric circles implying that vector qj can- 
not have any special significance as far as the DOS is 
concerned. As we shall see below the peaks at ±q7 arise 
solely due to the coherence factors, i.e. factors appearing 
in the numerator of Eqs. (01}. 

Integrals of the type appearing in Eq. (@J are familiar 
from the theory of massless relativistic Dirac fermions 
[2l| and can be evaluated most conveniently by exploiting 
the Feynman parameterization [2^. The exact result is 
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(5) 



For a given fixed energy uj function T[z) implies an in- 
verse square root singularity in both the real and imagi- 
nary parts of Aii n (q, uj) along an elliptic contour of con- 
stant energy given by E^ = 2uo (here E^ = 2q = 




a) non-magnetic 



b) magnetic 



FIG. 2: Interference patterns in the single-node approxima- 
tion, a) non-magnetic scattering Eq. 0, b) magnetic scat- 
tering Eq. (0. In both cases we use vf/va = 4. Note that 
the axes here represent the nodal coordinates qi (horizontal) 
and q2 (vertical), rotated by 45° relative to Fig. Q 



weighted by an angular factor {^ijq) 2 = (vaQ2/E cl ) 2 pro- 
ducing the largest amplitude at the two ends of the el- 
lipse, as illustrated in Fig. [3 These points of largest 
intensity coincide with ±q_7. We conclude that in this 
case the octet model works because of the special BCS 
coherence factors and not because of the DOS arguments. 
One can perform similar analyses for the internodal scat- 
tering and find a large response near some of the vectors 
q^ indicated in Fig. [I] Interestingly we find that the peak 
expected from the octet model near q4, which is absent 
in the experimental data [l5| , corresponds to an endpoint 
of a line of higher intensity. 

One may now inquire what is the physics leading to 
the angular modulation of the singularity displayed in 
Eq. and the consequent formation of the interference 
peaks in A(q, uj). A little thought reveals that the prin- 
cipal cause of this behavior is the particle-hole mixing. 
Charge is not a good quantum number in a superconduc- 
tor; non- magnetic impurities couple to charge and this 
leads to strong angular modulation of the quasiparticle 
response as one travels along the contour of constant en- 
ergy around the nodal point. To see this more clearly 
consider scattering by magnetic impurities. Magnetic 
impurities couple to quasiparticle spin and the latter re- 
mains a good quantum number in a superconductor. To 
describe this situation we replace T3 in Eq. J2J) by To = 1. 
This reverses the sign in front of the A + A_ term in Eq. 
(|3J). Such a sign reversal has profound consequences for 
the response function A(q, uj) which can be again evalu- 
ated in the low energy limit near a single node. In this 
case the result is 
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with A ~ A a high energy cutoff. The inverse square 
root singularity of Eq. (0 is replaced by a cusp along 
E^ = 2uj and, crucially, there is now no angular modu- 



\/v 2 F q 2 + v\q 2 ). More importantly this singularity is lation. The interference patterns for magnetic disorder 
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FIG. 3: Interference patterns in a DDW model of the pseu- 
dogap (top) compared with a dSC (bottom). In all figures 
uj 0.3A and A = D = 0.5t. 



will exhibit continuous lines near the center of the BZ 
(Fig. |2b) and will be qualitatively different from those 
with scalar disorder, as already noted in Ref. [T^ . 

The analysis of the low-energy limit of A(q, uj) thus 
shows that while the quasiparticle dispersion deter- 
mines the possible loci of strong interference, it is the 
coherence factors that determine the strength and select 
the actual location of the peaks on these loci. 

To conclude our discussion of the superconducting 
state we evaluate £n(q, uj) given by Eq. numerically 
with the full dispersion ek and parameters relevant to 
Bi-2212. In Fig. ^ we use the full single-impurity, inter- 
mediate phase shift T-matrix, T(ou) = To(uj)to-\-Ts(uj)ts^ 
as described in Ref. [2^ . We find that experimental pat- 
terns are best reproduced by a particular value of the 
phase shift that yields T%{uj) « —To (a;)*. This choice 
effectively kills the analog of the A + A_ term in Eq. J3J) 
and we are led to believe that in real systems this cancel- 
lation has a more fundamental cause, perhaps related to 
spatial fluctuations of the phase of A in a material with 
nanoscale electronic inhomogeneity 0, ^| . 

We now turn to the pseudogap state. If the pseudo- 
gap is due to an ordering in the p-h channel then the 
coherence factors will generally differ from those describ- 
ing the superconductor. It is clear from the previous 
discussion that such coherence factors will produce qual- 
itatively different interference patterns even if the DOS 
remains similar to that of a dSC. We now illustrate this 
general statement on the example of a d-density wave 
(DDW) state 12] that has been proposed to describe the 
pseudogap phase in cuprates. This is perhaps the most 
relevant example since for specially chosen parameters 
(/i = t' = 0) the DOS in DDW state is identical to that 
in dSC. We show below that even in this case FT-STS 
patterns are qualitatively different. 

Under the assumptions leading to Eq. one can show 
that for DDW state the interference pattern is given by 




FIG. 4: a) Interference patterns within the QED3 scenario 
with the anomalous dimension exponent r\ — 0.0, 0.2, 0.4 (left 
to right) and all other parameters as in Fig. [I] b) The FT-STS 
pattern calculated for a normal metal with the dispersion as 
in the dSC case and uj = 0.02£ (left) and uj = 0.045t (right). 



where the prime denotes summation over the reduced 
Brillouin zone. The DDW propagator reads Go(k, iuj) = 
[(iuj - e(J - e^r 3 - £>kT 2 ] -1 , with 4 = ^(e k + e k+Q ), 
d 7>(ek - e k +Q), Q (7r,7r) and DDW gap Ac 
(Cos k x — cos k y ). One obtains 
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D\ and c ±' de- 



with iQ± — iuj — e±, E± = y (e±) _ ~r ±^± 

noting k ± q/2 as before. At low energies the spectrum 
of a DDW quasiparticle is Dirac-like with the node fixed 
at (7r/2,7r/2). In this limit JSJ can be again evaluated 
analytically by performing the nodal approximation. We 
find that the result is given by Eq. © with uj replaced by 
uj + /i and va by vd, the slope of the DDW gap. Thus, 
we find that at low energies DDW state will produce 
FT-STS patterns similar to those expected for dSC with 
magnetic scattering, characterized by continuous lines in 
the Brillouin zone as opposed to the sharp peaks. This 
is confirmed by a full numerical evaluation of Eq. (|HJ) 
displayed in Fig. 03 The DDW patterns are markedly 
different from dSC even at half filling. 

We now consider the class of theories which describe 
the pseudogap as an incipient superconducting gap. In 
these theories the long range superconducting order is 
destroyed above T c but short range pairing correlations 
persist for T c < T < T*. Based on the above discussion 
we expect that such a pseudogap phase would inherit the 
FT-STS patterns characteristic of the superconducting 
phase since the coherence factors will locally retain their 
p-h character. As an example we consider the QED3 
theory of the pseudogap phase which has been proposed 
to describe a phase disordered dSC Q. In this theory 
fluctuating phase of the SC order parameter produces an 
emergent massless U(l) gauge field which mediates long 
range interactions between the fermions. As a result the 
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fermion propagator becomes incoherent |8|: 

G (k,zc)= [w2 + ^ + A 2 ]1 _, /2 . 0) 

and exhibits a Luttinger liquid like dynamics at long 
distances with the sharp quasiparticle poles replaced by 
branch cut singularities characterized by the anomalous 
dimension exponent n. The latter is believed to be a 
small positive number but its exact value is a matter of 
debate 8] . Here we treat rj as a parameter and show that 
the structure of FT-STS patterns is insensitive to its ex- 
act value. The QED3 propagator lacks the off-diagonal 
part, reflecting the absence of true SC long range order. 

One can again evaluate Eq. J2J) with the QED3 propa- 
gator ® analytically within the single node approxima- 
tion. For scalar disorder one finds a result similar to Eq. 
(0 with the square root singularity at = 2co replaced 
by a weaker l/(z 2 — l)^ _r? singularity. Thus, for n < | we 
expect patterns similar to those in dSC. The numerical 
results, presented in Fig. 21 indeed show that the QED3 
interference patterns retain peaks at the same positions 
as in the SC phase but the peaks become smeared for 
larger rj reflecting the incoherent nature of the fermionic 
excitations described by Eq. Q. 

Finally, Fig. H displays the FT-STS patterns for the 
normal metallic state (Ao = Dq = 0), which are expected 
to describe strongly overdoped cuprates. 

In conclusion, we have shown here that the quasi- 
particle interference patterns seen in the FT-STS re- 
veal the signatures of both the quasiparticle dispersion 
and, through their sensitivity to the quasiparticle coher- 
ence factors, the nature of the underlying electronic order 
present in the system. In particular, the latter deter- 
mines the basic characteristic features of the interference 



patterns. We have demonstrated, by general arguments 
and detailed calculations within several relevant models, 
that the superconducting order is very special in that it 
alone produces patterns consistent with the experimen- 
tal data. Based on this insight we have proposed a test 
for the nature of the pseudogap phase in cuprates using 
FT-STS. If the pseudogap is due to fluctuating SC or- 
der, then the FT-STS patterns above T c should remain 
qualitatively the same as those below T c . If, on the other 
hand, the pseudogap is due to static or fluctuating order 
in the p-h channel, such as SDW, CDW or DDW, the 
patterns above T c should be qualitatively different, gen- 
erally exhibiting continuous lines or peaks at different 
positions. 

As mentioned above preliminary experimental data 
on Bi-2212 [HI show FT-STS patterns above T c that 
are very similar to those found deep in the SC phase 
[IH ITH suggesting that the pseudogap state is of 
predominantly superconducting origin. One may ask how 
robust is this identification of the pseudogap state, pro- 
vided that the experimental data [X3 can be reproduced. 
Our analysis indicates that the p-p nature of supercon- 
ducting correlations plays critical role in formation of the 
patterns observed experimentally. Furthermore, despite 
significant effort, we were unable to construct a model 
with instability in p-h channel that would mimic these 
patterns. Therefore, we must conclude that data of Yaz- 
dani and co-workers |l7|, if correct, place very strong con- 
straints on the nature of the pseudogap state in cuprates. 
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Hoffman, D.-H. Lee, S. Sachdev, D.E. Sheehy, O. Vafek, 
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